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Abstract —This study aims to assess for a Stirling engine the influence of the overall heat transfer coeffi¬ 
cient variation on the optimum state and on the optimum distribution of the heat transfer surface con¬ 
ductance or area among the machine heat exchangers. The analysis is based on a Stirling machine 
optimization method, previously elaborated, which is now applied to a cycle with total heat regener¬ 
ation. The method was conceived for an irreversible cycle with heat transfer across temperature differ¬ 
ences at the source and the sink, and heat losses between the hot-end and the cold-end of the engine. 
Source and sink of finite thermal capacity as well as thermostats are considered. The new approach 
considers a linear variation of the overall heat transfer coefficient of the machine heat exchangers with 
respect to the local temperature difference. A comparison of the optimum state and the optimum distri¬ 
bution of the heat transfer surface conductance or area among the heater and the cooler is made for 
several cases. © 1998 Elsevier Science Ltd. All rights reserved 

Stirling engine Optimization method Irreversible cycle Heat transfer surface Conductance 
or area distribution 


NOMENCLATURE 


A = Heat transfer surface area 
C= Capacity flowrate 
h = Overall heat transfer coefficient 
hA = Heat transfer surface conductance 
k' = Proportional factor 
NTU = Number of heat transfer units 
Q = Heat transfer rate 
T= Temperature 
IU= Power delivery 
Greek characters 

<j> = Heat transfer rate at the source 
y = Ratio of specific heats ( c p /c v ) 

)/ = Cycle efficiency 
0 = Temperature ratio 
Subscripts 

c = Cold 

/= Working fluid of the cycle 
7/=Hot end, source 
i = Intermediary 
L = Cold end, sink 
loss = Related to heat losses 
min = Minimum 
r = Reduced 

S= Generally for the source or sink (as in T s ) 

T= Related to the total heat transfer surface conductance or area 
w = Warm 
0 = Reference 

1 = In 

2 = Out 


*To whom all correspondence should be addressed. Fax: 40-1-410.44.88; E-mail: mcos@theta.termo.pub.ro. 
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INTRODUCTION 

In the last few years, the Stirling cycle machine has been intensively studied in connection to its 
possible solar, extra-terrestrial and military applications. Recent papers have made major contri¬ 
butions to the study and optimization of the Stirling machine [1-4], 

Following the current trend, the authors have elaborated a Stirling engine optimization 
method [5], according to which the cycle was optimized with respect to the maximum power 
delivery. The method was conceived to provide the engine optimum state defined by the outlet 
temperatures of the sources, the temperatures of the working fluid during the isothermal heat 
input and output of the cycle, and the cycle efficiency. The optimum distribution of the heat 
transfer surface conductance or area among the engine heat exchangers were also obtained. 

The originality of this model came from the fact that the heat transfer rate available at the 
source is treated as a parameter, as opposed to past studies where the parameter was the source 
temperature. 

The authors extend here previous works on Stirling engines [6,7] by considering the influence 
of the overall heat transfer coefficient variation on the engine characteristics, particularly on the 
distribution of the heat transfer surface conductance or area among the heater and the cooler. 

MODELING AND OPTIMIZATION OF THE STIRLING ENGINE IRREVERSIBLE CYCLE 

The isothermal Stirling engine cycle with total heat regeneration is presented in Fig. 1. The 
continuous variation of the temperature difference between the source/sink and the cycle work¬ 
ing fluid suggests the finite thermal capacity of the source and the sink. The thermal reservoirs 
become thermostats when this AT —> oo. Also this case is analyzed in the paper. The external 
irreversibilities of the cycle are represented in this study by: 

—the heat input rate Qh across AT H = AT Hl — T w at the source and the heat output rate Q L 
across A T L = T Li — T c at the sink; 

—the heat loss, Q\ oss , between the hot-end (T Hl ) and the cold-end ( T Ll ) of the engine. 

The model developed by Costea [6] is extended here to show the influence of the overall heat 
transfer coefficient variation on the performance of the engine. Generally, the overall heat trans¬ 
fer coefficient of the heat exchanger is not constant with respect to the temperature variation. 
Actually, the overall heat transfer coefficient may be considered as any function of AT due to 
the isothermal heat transfer input/output in a Stirling machine. 

The general case of the non-linear dependency of the overall heat transfer coefficient, h, on 
the temperature variation is difficult to model, despite the fact that any variation can be brought 



Fig. 1. Temperature-entropy diagram of a Stirling engine cycle with source and sink of finite thermal 

capacity. 
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Table 1. Presentation of the studied cases in the present optimization of the Stirling engine 


extf corresponds to C ex t = min(C ex t>C), where C ext = (Ch,Cl), with the assumption Cjj= min(C//,C/.) 
intj corresponds to C = min(C cx t,C), and a, [I, S are given by the following constraints: 


(a) Fixed heat transfer surface conductance of the heater and the cooler 

I’loAh, IiloAl = fixed. 

(/I) Fixed total heat transfer surface conductance 

h m A H + h L0 A L = fixed = h n A T . 

(<5) Fixed total heat transfer surface area 


to a linear one within small intervals. Colburn, for example, assumed a linear variation of h 
with respect to the local temperature difference in the heat exchanger [8], as follows: 

hsi = hsio • (1 + h'si ■ kT St ) (1) 

where h S m is a reference value for the overall heat transfer coefficient of the source/sink heat 
exchanger, and h' Si is the proportional factor of this variation. 

The authors of the present study adopted the above assumption as a new approach in Stirling 
engine theoretical investigation. 

The heat losses between the source and the sink are included in the analysis via an additional 
heat transfer conductance, and is expressed by equation (2): 

gloss = (hA\ oss .(T m - T u ). (2) 

The mathematical model is based on the First and Second Law of Thermodynamics applied to 
the cycle, the energy balance at the source, and heat transfer equations. All the aforementioned 
irreversibilities are included. The model consists of a nonlinear system of equations (objective 
function and two constrains). Solutions to the model are found by using Lagrange’s multiplier 
method. 

For each set of assumptions (see Table 1), the model yields different expressions for the 
objective functions and the constraints. Two representative cases were chosen to illustrate the 
resulting system of equations: 

1. Source and the sink of finite thermal capacity: FTC flux WL-ext a (see Table 1). 

(i) The objective function: 


y Ch ao [ 1 _ exp(—NTUj/p) 1 / 1 - A0 U \ 

C 1 + k' [fr A() m [ 1 - exp(—NTUtf 0 )]J V 6 m - A0 m )' 


( 3 ) 


(ii) The constraints: 

• the cycle entropy generation 

2£h\Q 1 _ exp(—NTU//o) _ AT Ad 1_ exp(—NTU m) _ 

C l+F„A0 H1 [l-exp(-NTU„o)]J C l+F^A^fl-expf-NTUio)]] 

Om - A d m + 1 - A 0 Ll 


( 4 ) 
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the energy balance at the source 


1C H 

c 


AB 


HI 


1 - 


exp(— NTU//o) 


l+k' Hr AB m [\ - exp(-NTU^o)] 


4>r — NTU loss (0//i - 1) 


= 0. (5) 


(2) Source and the sink of constant temperature: CT-temp-ext-a. 
(i) The objective function: 


A B, 


W, = -I^ntUhoO +VhA<)h)Mh(i -2 _ Afe 

(ii) The constraint: 

• the cycle entropy generation 

^NTUM,A0tf (1 + k' Hr A0 H ) ^NTU lo A 0 l (1 + k' Lr AB L ) 


( 6 ) 


Bh-AB 


■ + ■ 


H 


1 - ABj 


= 0 . 


(7) 


Note that both systems are representative for the case with fixed heat transfer surface conduc¬ 
tance corresponding to the reference value of the overall heat transfer coefficient (7i m A H and 
Jilo^l) °f the heater and the cooler. For each of the other two cases—and S —the NTU// 0 
and the NTU^o will be modified corresponding to the new conditions imposed by the new 
constraint. Also a new variable (heat transfer surface conductance or area) will be added to 
those of (a)-case. 

The non-dimensional variables of the model are: 

—the temperatures: 


FTC A B Ll = 1 - AB m = Tm rr Tw 

1 L\ 


Bm = 


’Ll 


T m 


’ L\ 


-the heat transfer surface conductances or areas: 

hm^n BlqAl 


08) 


BtoAt 


hroAr 


® % 


Al 

At 


The model parameters are: 

—the total reference heat transfer conductance equivalent NTU: 

a , . TTT , BtqAt . m-ttt hroAr 

<>ext : NTU 7 X) — —7 <>mt : NTU™ — — 7 — 

Cc C 

—the proportional factor of the overall heat transfer coefficient variation: 


h 




7 1 77 

K = h»' Tu 


-ratio of the reference overall heat transfer coefficient of the sink and the source: 

(<5) rh = — 
hm 


-the heat loss conductance equivalent NTU: 

(^)loss 


<>ext : NTUi oss = 


C c 


Oint : NTUi oss = 


( hA ) 


loss 


c 


-the heat transfer rate available to the source: 

FTC ([>,. = 




C/y • T li 


= fir. 
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Fig. 2. Influence of the heat transfer rate available at the source on the Stirling engine optimum tem¬ 
peratures for finite thermal capacity and constant temperature of the source 

(9 h = 3; NTUtd = 5; NTU loss = 0.5; h Hr = 0.1; J,, = 0.1). 


Note that T [A from the above FTC case equations becomes T L in those corresponding to the 
CT case. 


RESULTS AND DISCUSSION 

The summary of the possible cases to be studied with this optimization model is presented in 
Table 1. Analytical results were obtained only for the cases without heat losses between the 
source and the sink, and they show that there is no physical optimum. 

Note that when the source temperature is fixed, the energy balance at the source cannot be 
written anymore. Since it was the only equation to account on the heat losses between the 
source and the sink, there is no more meaning to distinguish two cases corresponding to the 
heat losses. Just to complete the optimization results, one can calculate separately the total heat 
transfer rate (Q H + Qi oss ) supplied from the source. 

Also there is no meaning to consider C = min(C, C ex t) for the case of finite thermal capacity 
of the source and the sink because of the isothermal processes of cycle heat input/output. 
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Fig. 3. Influence of the heat transfer rate available at the source on the Stirling engine power output, 
cycle efficiency, and optimum heat transfer surface conductance for finite thermal capacity and constant 
temperature of the source 

(0 H = 3; NTUnj = 5; NTU loss = 0.5; h Hr = 0.1; h Lr = 0.1). 


A parametric study was conducted in order to observe the effect of the overall heat transfer 
coefficient variation on the optimum state, and on the optimum distribution of the heat transfer 
surface conductance or area among the engine heat exchangers. 

We attempted to provide realistic central values for the model parameters, as well as for their 
variation range. In order to compare the results with previous data obtained for Stirling cycle 
solar application, the central value of the heat transfer rate available to the source was chosen 
high. The present investigation considered the following parameter ranges: 

— 4> r = ( 0.1, 10), central value 0.75 or 0 m (1.5, 8), central value 3; 

—NTUjp = (0.1, 5), central value 3; 

—NTUi oss = (0, 5), central value 0.5; 

— h'Hr, h' Lr = (0 , 1), central value 0.1; 

—rh = ( 0.01, 4), central value 0.1. 

Some selected results are graphically presented in Figs 2-5. 

Figures 2 and 3 show that, for the case of source and sink of constant temperature, the tem¬ 
perature differences, A9* m , A9* L \, decrease and thus the power output will increase relative to 
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Ratio of reduced overall heat 
transfert coefficient of the sink 
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Ratio of the reduced overall 
heat transfer coefficient of the 
sink and the source 


Fig. 4. Comparison of the engine temperatures for the case of finite thermal capacity of the source and 
sink, with heat transfer rate available at the source or fixed inlet source temperature 

(fir = 0.75 or 0m = 3; NTUtd = 5; NTUi oss = 0.5; h Hr = 0.1; rh = 0.1). 


the corresponding values of the FTC case. The cycle efficiency and the working fluid cold-end 
temperature are not affected. Significant differences between the results of these cases are regis¬ 
tered for high value of 4> r . However, for small values of <j) r , only the heat transfer conductances 
of the heater and cooler are significantly sensitive with respect to the heat transfer rate variation 
(decreasing from 2.5 to 1.7 for the heater, increasing from 1.5 to 2.65 for the cooler, when r/),. 
increase from 0.1 to 1). 

Results in Figs 4 and 5 show that an imposed heat transfer rate at the source yields smaller 
value of the source inlet temperature comparatively to 9 H = 3 of the FTC-temp case. As a con¬ 
sequence, the power output of the engine decreases drastically by about 300%. The associated 
heat transfer surface area distribution was expected, since an increase in rh means h L0 f and 
JihoU namely, the opposite variation of the corresponding areas. 

Remarkably, there is no difference between the heat transfer surface area distribution for 
these two cases. As a final remark on Fig. 3, it is obvious that is worthless to enhance very 
much the overall heat transfer coefficient of the cooler, as the important variation of the cycle 
characteristics are registered for r/7 = (0.01, 2). 
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Ratio of reduced overall heat 
transfert coefficient of the sink 
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Fig. 5. Comparison of the engine power output, cycle efficiency, and optimum heat transfer surface 
area distribution for the case of finite thermal capacity of the source and sink, with heat transfer rate 
available at the source or fixed inlet source temperature 

(fir = 0.75 or d H \ = 3; NTUn, = 5; NTU ioss = 0.5; h Hr = 0.1; rh = 0.1). 


CONCLUSIONS 

A theoretical model was presented for the Stirling engine cycle optimization based on a linear 
variation of the overall heat transfer coefficient with the temperature difference. Solutions corre¬ 
sponding to different sets of assumptions and parameter values were used for a comparison of 
the engine characteristics. Results pointed out either optimal variation range for some model 
parameters, or some significant differences of the power output, source and sink temperature 
difference, heat transfer characteristics values with respect to each of the studied cases. The 
model allows further investigation that can provide additional insight into the optimal design of 
a Stirling engine, being currently developed by the authors. 
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